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Abstract 

The dynamics of the Einstein- Vlasov equations for a class of cosmological models 
with four Killing vectors is discussed in the case of massive particles. It is shown that 
in all models analysed the solutions with massive particles are asymptotic to solutions 
with massless particles at early times. It is also shown that in Bianchi types I and II 
the solutions with massive particles are asymptotic to dust solutions at late times. That 
Bianchi type III models are also asymptotic to dust solutions at late times is consistent 
with our results but is not established by them. 



1 Introduction 

The most popular matter content by far in the study of spatially homogeneous cosmological 
models is a perfect fluid with linear equation of state (see e.g., the book 12]). It is important 



to know if the results obtained for this class are structurally stable if we change the matter 
content. Thus it is of interest to investigate other types of sources. Here we will consider 
certain diagonalizable locally rotationally symmetric (LRS) spatially homogeneous models 
with collisionless matter. This class of models was previously studied in the case of massless 
particles in |10[ |. Here we will focus on the case with massive particles. We will recast 
Einstein's field equations into a form so that one part of the boundary of the state space for 
the massive case can be identified with the state space for the massless case while another 
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part can be identified with the state space for the corresponding dust equations. (In addition 
other parts of the boundary have the interpretation of state spaces associated with certain 
models with distributional matter.) It will be shown that these boundary submanifolds are 
intimately connected with the early and late time behaviour of the LRS massive collisionless 
gas models respectively. 

The results of our analysis can be summarized as follows. Consideration is restricted to 
models of Bianchi types I, II and III. This is enough to display a large variety of phenomena. 
At early times, i.e. close to the initial singularity, the dynamics of solutions with massive 
particles mimics closely the dynamics for the corresponding symmetry type with massless 
particles. In particular there are solutions whose behaviour near the singularity is quite 
different from that of any fluid model of any of these Bianchi types. At late times, i.e. in a 
phase of unlimited expansion, the general picture is that the dynamics resembles that of a dust 
model. This is proved for Bianchi types I and II. For type III the results are consistent with 
dust-like asymptotics but we were not able to prove that this is what happens. If kinetic 
theory with massive particles always behaved like dust at late times this would provide a 
justification of the use of a fluid model in that regime. 

The outline of the paper is as follows. In section 2 we derive the dynamical system. 
Sections 3, 4 and 5 analyse the models of types I, II and III respectively, with the main 
results being stated in Theorems 2.1, 3.1 and 4.1. In section 6 we conclude with some 
remarks and speculations. An appendix contains some information about dynamical systems 
which is applied frequently in the paper. 

2 A dynamical systems formulation 

We will consider LRS models for which the metric can be written in the form 

ds 2 = -dt 2 + gnitW 1 ) 2 + g 2 2(t)((9 2 ) 2 + (fl 3 ) 2 ) , (1) 

where Q % are suitable one- forms describing the various symmetry types. The energy- momentum 
tensor Tij for the Einstein- Vlasov system with massive particles is assumed to be diagonal 
and is described by 

P = J fo(v i )(m 2 +g 1 \v 1 ) 2 +g 22 ((v 2 ) 2 + (v 3 ) 2 )) 1 / 2 (detg)- 1 / 2 dv 1 dv 2 dv 3 , 
Pi = J h(vi)g U (vi) 2 (m 2 + g l \ Vl ) 2 + g 22 ((v 2 ) 2 + (v^y^idetgy^dvrdvidvz , (2) 

where p is the energy density and pi = T % i the pressure components of the energy-momentum 
tensor. The function fo is determined at some fixed time to by fo(vi) = f{to,vi) where / is 
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the phase space density of particles. The covariant components Uj are independent of time. 
The function /q satisfies the condition fo(vi,V2,V3) = F(v\, (v 2 ) 2 + {v 3 ) 2 )- 

Some further technical conditions will be imposed on /q. It is assumed to be non-negative 
and have compact support. It is also assumed that the support does not intersect the co- 
ordinate planes Vi = 0. A function /o with this property will be said to have split support. 
The reason for the assumption of split support will be seen later. In the following it will 
always be assumed without further comment that the data considered have split support. 
It follows from the assumptions already made that fo(xi) = fo(-Xi) for i = 2,3. It will 
be assumed that this also holds for i = 1 and functions fo with this property will be called 
reflection-symmetric. This ensures that the form of the phase space density of particles is 
compatible with a diagonal metric and, in particular, that the energy-momentum tensor is 
diagonal. For the symmetry types to be considered in the following it then follows that 
the entire system consisting of geometry and matter is invariant under three commuting 
reflections. For this reason, solutions where the metric is diagonal and /o has the symme- 
try properties just mentioned will be called reflection-symmetric. A solution is said to be 
isotropic if fo(v 1 ,v 2 ,v 3 ) = F((^i) 2 + (v 2 ) 2 + (v 3 ) 2 ) and if g n oc g 22 oc g 33 for all time. 

The momentum constraints are automatically satisfied for these models. Only the Hamil- 
tonian constraint and the evolution equations are left. Instead of considering a set of second 
order equations in terms of e.g., a and b, where 

a 2 = gn , b 2 = g 22 , (3) 

we will reformulate these equations as a first order system of ODEs by introducing a new set 
of variables. The mean curvature tik (where kij is the second fundamental form) is given by 

txk = -{ar l da/dt + 2b~ 1 db/dt) . (4) 

A new dimensionless time coordinate r is defined by — | f t trk(t)dt for some arbitrary fixed 
time to- (We will follow the conventions in filf . The time variable thus differs by a factor 
3 from the one in |l0|]). In the following a dot over a quantity denotes its derivative with 
respect to r. The Hubble variable H is given by H = — tr/c/3. Now define the following 
dimensionless variables: 

m 2 j(cr 2 + 2£T 2 +m 2 ) , 
b 2 /(b 2 + 2a 2 ) , 
a 2 (a 2 /b 4 )(tvk)- 2 , 
3a 3 b~ 2 (trky 2 , 

-Zib^db/dtjitrk)- 1 - 1 , (5) 

where a 2 is 1 for Bianchi types II, VIII, IX and for Bianchi types I, III and the Kantowski- 
Sachs (KS) models. The coefficient a 3 is 1 for types III and VIII. It is —1 for KS and type 



z = 

s = 

M 2 = 

M 3 = 

S + = 
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IX and for types I, II These variables lead to a decoupling of the equation for the only 
remaining dimensional variable H (or equivalently trfc) 

H = -(l + q)H, (6) 

where the deceleration parameter q is given by 

q = 2T, 2 + + + R) . (7) 

The quantity R is defined by 



R=(pi + 2 P2 )/p , (8) 



where 



pi/p = (l-z)sg 1 /h, 

Pz/p = |(1 - z)(l - s)g 2 /h , 

9i,2 = I fo(vi)(vi, 2 ) 2 [z + (1 - z)(s( Vl ) 2 + |(1 - s)((?; 2 ) 2 + (« 3 ) 2 ))]- 1/2 d«id«2d« 3 , 



h = / / («i)[« + (1 - z)(s( Vl ) 2 + i(l - s)((v 2 ) 2 + ( W3 ) 2 ))] 1/2 ^i^ 2 ^3 • (9) 



The assumption of split support ensures that the function R(s,z) is a smooth (C°°) 
function of its arguments. The related quantity itl+ defined by 

R+ = (p 2 -Pi)/p ■ (10) 

is a smooth function of s and z for the same reason. 

The normalized energy density Q = p/(3H 2 ) is determined by the Hamiltonian constraint 
and, in units where G = l/8ir, is given by 

n = 1 - T, 2 + - M 2 - M 3 . (11) 

The assumption of a distribution of massive particles with non-negative mass leads to 
inequalities for R, R + and S7. Firstly, < R < 1 with R = only when z = 1 and R = 1 
only when z = 0. Secondly, — i? < i? + < ^ii with i? + = ^i? for s = and i?+ = — i? for 
s = 1. Thirdly > 0. Using these inequalities in equation (0) in turn results in < q < 2 
for Bianchi types I, II, III and VIII (i.e., the same inequality as for causal perfect fluids, see 
Pi). 



The motivation for the variable s comes from more general diagonal models where it is convenient to 
introduce variables of the type Sj = g t '/(g 11 + g 22 + g 33 )- s is simply si in the case when g 22 = g 33 . 
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The remaining dimensionless coupled system is: 



£+ = -(2 -g)S + -S + + QR+ 

s = 6s(l — , 

z = 2z(l-z)(l + E + -3E + 8) 

M 2 = 2(q - 4S+)M 2 , 

M 3 = 2(5 - £+)M 3 , 



where 5+ is given by 



5+ = -4M 2 - M 3 . 

There are a variety of submanifolds corresponding to different symmetry types: 



(12) 
(13) 



Si 


M 2 


= M 3 = , 








Sn 


M 2 


> , M 3 = 








Sin 


M 2 


= , M 3 > 








Sks 


M 2 


= , M 3 < 








Syin 


M 2 


> , M 3 > 


(1- 


s)M 3 


= 6M 2 s 


Six 


M 2 


> , M 3 < 


(1- 


s)M 3 


= -6M 2 s 



(14) 

The relationship between the various models can be visualized in a symmetry reduction 
diagram given in Fig. |l| (a collective treatment of the corresponding vacuum models from a 
Hamiltonian perspective and with the aim of quantizing the models was given in ||]). Note 
that while this diagram accurately reflects the relationship of the geometry in the different 
cases, the relationship of the matter content is more subtle when types VIII or IX are involved. 
This complication does not occur for the Bianchi types studied in detail in this paper and 
will therefore not be discussed further here. 

Note that a non-negative energy density implies that > 0, which in turn implies that our 
variables are bounded for types I, II, III and VIII, since M 2 and M 3 are non-negative and since 
by definition z and s are bounded. These models expand indefinitely. The KS and type IX 
models are recollapsing models and since H becomes zero at the point of maximal expansion, 
the Hubble- normalized variables blow up at this point. However, one can find other variables 
that are bounded along the lines found in fnj. Neither are the above variables 'optimal' 
for the other LRS models. One can adapt to the particular mathematical features these 
models exhibit. However, we choose to use the above formulation since the present variables 
are easier to interpret physically and are naturally generalizable to more general non-LRS 
models. For simplicity we will from now on study Bianchi types I, II, and III. 
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VIII IX 




III II KS 




Figure 1: Symmetry reduction diagram for the diagonal LRS models. 



It is of interest to note that the metric functions a, b are expressible in terms of s, z in the 
massive case. The relations are 

a 2 = z(m 2 s(l- z))-\ b 2 = 2z(m 2 (l-s)(l- z))- 1 . (15) 

In addition to the symmetry submanifolds there are a number of other boundary sub- 
manifolds: 

z = 0,1, 
s = 0, 1 , 

n = o . (16) 

The submanifold z = corresponds to the massless case. The submanifold z = 1 leads to a 
decoupling of the s-equation, leaving a system identical to the corresponding dust equations. 
The submanifolds s = 0, s = 1 correspond to problems with fo being a distribution while 
O = constitutes the vacuum submanifold with test matter. Apart from these solutions 
there exists an isotropic solution in Bianchi type I characterized by S + = R + = and a 
constant value for s that depends on the function /q. 

Including these boundaries yields compact state spaces for types 1,11 and III. In order to 
apply the standard theory of dynamical systems the coefficients must be C 1 on the entire 
compact state space G of a given model. This is necessary even for uniqueness. In the 
present case it suffices to show that R and -R+ are C 1 on G, i.e, that they are C 1 for s,z 
when 0<s<l,0<z<l. As has already been pointed out, this follows from the assumption 
of split support, which even implies the analogous statement with C 1 replaced by C°°. It 
would be possible to get C 1 regularity under the weaker assumption that /o vanishes as fast 
as a sufficiently high power of the distance to the coordinate planes. We have not, however, 
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examined in detail how large the power would have to be since this is of little relevance to 
our main concerns in this paper. 

Of key importance is the existence of a monotone function in the 'massive' interior part 
of the state space: 

M = (s(l- s) 2 )- 1/3 z(l- z)- 1 , 

M = 2M . (17) 

Note that the volume ab 2 is proportional to M 3 / 2 . This monotone function rules out any 
interior uj- and a- limit sets and forces these sets to lie on the s = 0, s = 1, 2 = or 2 = 1 
parts of the boundary. 

3 Type I models 

It is natural to start investigating the type I system since it is a submanifold of the state 
space of all other symmetry types. The physical state space, G, of these models is given by 
the region in R 3 defined by the inequalities — 1 < S + < 1, < s < 1 and < z < 1. 

To understand the dynamics of the type I models, it is necessary to determine the station- 
ary points and their stability. The coordinates, in terms of (£+, s, z), of the various stationary 
points are the following: (0, «„, 0),(±, 0, 0), (1, 0, 0),(1, 1, 0),(-l, 1, 0),(1, 0, 1),(1, 1, 1),(-1, 1, 1), 
where sq is a particular constant value of s depending on the function / (see |l(]]). These 
points are called Pi, ...P$- (Note that they are numbered differently in the massless case com- 
pared to those in Jl0| .) In addition there exist two lines of equilibrium points, (—1, 0, K), (0, F, 1), 
denoted by Li,L2, where K and F are constant values. The points Pi, P2, P4, P&-, P7, Ps are 
hyperbolic saddles while P5 is degenerate, with one zero eigenvalue. The point P3 is a 
hyperbolic source. The line L\ is a transversally hyperbolic saddle while the line L2 is a 
transversally hyperbolic sink. (For an explanation of this terminology we refer to the ap- 
pendix.) 

The state space together with equilibrium points and separatrix orbits is depicted in Fig. 

& 

The main result in this section is the following theorem: 
Theorem 3.1 If a smooth non- vacuum reflection-symmetric LRS solution of Bianchi type 
I of the Einstein- Vlasov equations for massive particles is represented as a solution of ( |l2|) 
with M2 = M3 = then for r — > 00 it converges to a point of the line L2. For r — > — 00 there 
exists 

(i) a single (isotropic) solution that converges to Pi and 

(ii) a one-parameter set of solutions lying on the unstable manifold of P2 and 

(iii) all remaining solutions belong to a two-parameter set (the generic case) of solutions 
converging to P3. 
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Figure 2: The LRS type I state space together with equilibrium points and separatrix orbits. 

This will be proved in a series of lemmas. We refer to for terminology from the 

theory of dynamical systems. 

Lemma 3.1 There exist open neighbourhoods U\ and U% of the point P3 and the line L2 
respectively such that: 

(i) if a solution belongs to U\ at any time it belongs to U\ at all earlier times and its a-limit 
set consists of the point P3 alone. 

(ii) if a solution belongs to U2 at any time it belongs to U2 at all later times and its w-limit 
point consists of a single point of the line L2. 

Proof Part (i) follows from the fact that P3 is a hyperbolic source and the Hartman-Grobman 
theorem. Part (ii) follows from the fact that L2 is a transversally hyperbolic sink and the 
reduction theorem ( fLOU , Theorem Al). 

As a step towards analysing the dynamics of the full system we determine the w-limit 
points of solutions of the dynamical system on the parts of the boundary of G defined by 
s = and s = 1. This information will later be combined with the monotone function M 
when determining the w-limit sets of solutions of the full system. In the case of the a-limit 
sets the monotone function alone accomplishes the same thing. 

Lemma 3.2 A solution of the restriction of the system to the part of the boundary of G 
defined by s = 1 for which neither z nor S + take on one of their limiting values has the 
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endpoint of L2 with s = 1 as its u-limit set. 

Proof If S + > at any time, then S + is decreasing. The rate of decrease remains uniform 
as long as S + does not tend to zero. It follows that after a finite time S + must be strictly 
less than 1/2. On the other hand, z is monotone increasing in the region S + < 1/2 and the 
rate of increase remains uniform as long as z does not tend to one. It follows that z — > 1 as 
r — > 00. If E_|_ tends to zero in this limit then the conclusion of the lemma holds. Otherwise 
S + must become negative at some time. Thus it can be seen that any w-limit points satisfy 
z = 1 and —1 < S + < 0. From part (ii) of Lemma 3.1 it follows that any solution which 
enters U2 has the desired w-limit set. Since the w-limit set is a union of orbits, it is possible 
as a consequence to exclude the points with z = 1 and — 1 < S + < from the w-limit set. To 
complete the proof of the lemma it remains only to exclude the point Pg from the cj-limit set. 
This point is a hyperbolic saddle of the restriction of the system to s = 1 and so it follows 
from the discussion in the appendix and what has been proved already that it cannot belong 
to the w-limit set. For if P$ belonged to the w-limit set points of its stable and unstable 
manifolds would also have to do so, and this has already been ruled out. 
Lemma 3.3 A solution of the restriction of the system to the part of the boundary of G 
defined by s = for which neither z nor S + take on one of their limiting values has the 
endpoint of L2 with s = as its w-limit set. 

Proof Along any solution of this system z is monotone increasing on the part of the state 
space of the restricted dynamical system with S + 7^ —1 and z(l — z) 7^ 0. Hence, by the 
monotonicity principle (see [13|), any cj-limit point must satisfy z = 1 or S + = — 1. However 
S + is increasing for S + close to but not equal to —1. Hence there can be no cj-limit points 
with = — 1. It follows that z tends to one as r — > 00 for any solution and any w-limit 
point satisfies z = 1. From Lemma 3.1, any solution which enters XJ% has the desired w-limit 
set. Arguing as in the proof of Lemma 3.3 allows points with —1 < S + < and < S + < 1 
to be excluded. The point Pq, which is a hyperbolic saddle of the restricted system, can be 
eliminated in the same way as was done in the case of Pg in the proof of Lemma 3.2 using the 
results of the discussion in the appendix. Finally, the non-existence of cj-limit points with 
S + = — 1, already mentioned above, shows that the endpoint of the line L\ cannot lie in the 
w-limit set. 

Lemma 3.4 If a solution lies in the interior of G, then unless it lies on the unstable manifold 
of Pi or P2 its a-limit set consists of the point P3. 

Proof Consider a solution in the interior of G which does not lie on the unstable manifold 
of Pi or Pj. If it intersects U\ then by Lemma 3.1 its a-limit set consists of the point P3. 
There can be no other a-limit points in U\. Because the function M tends to zero along the 
solution as t — > — 00 the a-limit set must be contained in the surface z = 0. Recall that the 
surface 2 = corresponds to the case of massless particles which was analysed completely in 



10[. (Note that the stationary points were numbered differently in that paper.) Consider 



the boundary of the surface z = 0. Arguing as in the proof of Lemma 3.2, the lines joining 



9 



i-*3 to Pi and P4 can be excluded from the a-limit set. The discussion of the appendix and 
the fact that P4 is a hyperbolic saddle with stable manifold S + = 1 and unstable manifold 
the line connecting P4 to P5 can be used to exclude that line and the point P4 itself. The 
line connecting P5 to the endpoint of the line L\ can be excluded in an analogous way, noting 
that the non-hyperbolic point P5 is also covered by the discussion of the appendix. The point 
P5 is also excluded by this argument. Applying the reduction theorem allows the line joining 
the endpoint of the line L\ to Pi to be excluded together with the endpoint of L\. At this 
stage we can also exclude the point Pi itself, using the results of the appendix again and the 
fact that by assumption the solution does not lie on the unstable manifold of Pi. Thus the 
only point of the boundary of the set z = which can belong to the a-limit set is P3. Now 
suppose that a point of the interior of the surface belongs to the a-limit set. If it is a point 
of the unstable manifold of Pi then Pi also belongs to the a-limit set, in contradiction to 
what has just been proved. If it is some other point other than Pi then, using the fact that 
the a-limit set is a union of orbits and Theorem 3.1 of [pTfl , it follows that P3 belongs to the 
a-limit set and we obtain a contradiction again. Finally, if it were Pi then the results of the 
appendix would imply that other points of the interior would belong to the a-limit set, and 
this has just been ruled out. 

Lemma 3.5 The w-limit point of each solution in the interior of G is a point of the line Li. 
Proof Note first that the function M goes to infinity along any such solution as r — > 00. 
It follows that any w-limit point must satisfy z = 1, s = or s = 1. If the solution enters 
the set Ui then by part (ii) of Lemma 3.1 the w-limit set is as claimed. There are no other 
w-limit points of any solution in 17%. Consider now the evolution of £+ on the surface z = 1. 
It either increases from —1 to or decreases from 1 to 0. Since the w-limit set is a union of 
orbits, we conclude that no point of the interior of the surface z = 1 or its boundary lines 
s = and s = 1 other than the points of the line Li can belong to the o;-limit set. Using once 
more the fact that the o;-limit set is a union of orbits, it is possible to exclude the interior 
of the surface s = 1 from the cj-limit set by Lemma 3.2 and the interior of s = by Lemma 
3.3. Now all remaining possibilities other than points on Li will be excluded successively. 
The nature of the line L\ as a transversely hyperbolic saddle suffices to eliminate it, as well 
as the lines joining it to Ps and Pi. The point P3, being a hyperbolic source, is clearly ruled 
out, and with it the lines joining it to Pi and P§. Further applications of the results of the 
appendix rule out the remaining lines, namely those joining P$ to P5, P5 to P4, P4 to P7 and 
P7 to Pq. It follows that the w-limit set is contained in the line Li. Applying the reduction 
theorem then shows that the w-limit set is a single point of Li. 

The results of Lemma 3.4 and Lemma 3.5 together imply Theorem 3.1. 

Theorem 3.1 has been formulated entirely in terms of the dynamical systems picture. It 
should, however, be pointed out that this allows asymptotic expansions for all quantities of 
geometrical or physical interest near the singularity or in an expanding phase to be obtained 
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if desired. For example, in an expanding phase in type I the following expansions can be 
derived: 





(18) 


s = s - |s (l - so)^ 1 + o{t~ l ) 


(19) 


Z = 1-^-4/3 + ^-4/3) 


(20) 


H = fi" 1 + 0(t~ 7 / 3 ) 


(21) 


n = 1 - at' 2 + o(i~ 2 ) 


(22) 


p = |t- 2 -|a 2 r 4 + o(t- 4 ) 


(23) 


Pl = o(t- 10 / 3 ) 


(24) 



Here a and (5 are constants depending on the solution. It should be emphasized that these are 
not just formal expansions, but rigorous results which emerge from the dynamical systems 
analysis. 

A particular consequence of Theorem 3.1 is that all LRS type I models isotropize at late 
times. This was already proved by other means in ||, where it was also shown that non-LRS 
models of Bianchi type I isotropize and have dust-like behaviour for r — > oo. 

4 Type II models 

The physical state space, G, of the LRS type II models is given by the region in R 4 defined 
by the inequalities M 2 > 0, < s < 1, < z < 1, and 1 - - M 2 > 0. 

The coordinates, in terms of (£+, s, z, M%), of the various stationary points are the follow- 
ing: (0, So , 0,0),(i, 0,0,0), (1, 0, 0, 0),(1, 1, 0, 0),(-l, 1, 0, 0),(1, 0, 1, 0),(1, 1, 1, 0),(-l, 1, 1, 0), 
(ijj, 1, 0, <§;),(|, 0, 1, j 1, 1, (g), where so is the same particular constant value of s that 

appeared in the previous type I section. These points are called Pi, ...Pn (note that they are 
numbered differently than in JlO[] , in the massless case). In addition there exist two lines of 
equilibrium points, (—1, 0, K, 0), (0, F, 1, 0), denoted by Li, L2, where K and F are constants. 
The first eight stationary points and the two lines correspond to points and lines of the same 
name in the Bianchi I system and their coordinates are obtained by appending a zero to those 
of the Bianchi I points. The points Pi, P 2 , P3, P4, Pg, P7, Pg, P9, Pro are hyperbolic saddles 
while P5 is degenerate, with one zero eigenvalue. The point Pn is a hyperbolic sink with two 
real and two complex eigenvalues. The lines L\ and L 2 are transversally hyperbolic saddles. 

To prove results about the global properties of solutions it is helpful to use certain mono- 
tone functions. The first is defined for s < 1 by 

Z x = (2s/(l - s)) 4 / 3 M 2 (25) 
11 



This is obtained by rewriting the function whose time derivative was calculated in equation 
(23) of [10] in terms of the variables of this paper and observing that it remains monotone 
in the massive case. It satisfies Z\ = 2qZ\. The second is obtained by combining Z\ with 
the monotone function M available for all the Bianchi types considered in this paper. Let 
Z 2 = ZiM~ 2 = 2 4 / 3 s 2 M 2 z- 2 (l - z) 2 for z > 0. It satisfies Z 2 = 2(q - 2)Z 2 . The function Z x 
is defined on the part of the Bianchi II state space where s ^ 1 and monotonically increasing 
except where it vanishes. This is clear if q ^ 0. If q = it follows that S + = and M 2 = 1 
and at points satisfying these conditions X+ ^ 0. The function Z 2 is defined on the part 
of the Bianchi II state space where z > and is monotonically decreasing except on the set 
where it vanishes, since q = 2 implies Z 2 = 0. 

Theorem 4.1 If a smooth non- vacuum reflection-symmetric LRS solution of Bianchi type 
II of the Einstein- Vlasov equations for massive particles is represented as a solution of (12) 
with M3 = 0, then for r — > 00 it converges to P\\. For r — > —00 there exists 

(i) a one-parameter set of solutions converging to the unstable manifold of Pi and 

(ii) a three-parameter set of all remaining solutions converging to the heteroclinic cycle on 
the z = submanifold, consisting of the orbits connecting the z = endpoint of the line L\ 
to P5, P5 to P4, P4 to P3 on the type I boundary and P3 to the z = endpoint of the line L\ 
via the vacuum boundary. 

Lemma 4.1 If a solution belongs to the interior of the type II state space then any a-limit 
point satisfies z = and sM 2 = 0. Any w-limit point satisfies s = 1 and {z — 1)M 2 = 0. 
Proof From the evolution equation for M it follows that z = for any a-limit point and that 
for any cj-limit point z = 1, s = or s = 1. Next the monotonicity principle will be applied 
to the functions Z\ and Z 2 . Applying it to Z\ on the region where Z\ ^ shows that for 
any a-limit point s = or M 2 = 0. It also shows that there are no u;-limit points with s = 0. 
Combining this with the information obtained already shows that any cj-limit point satisfies 
z = 1 or s = 1. If z ^ 1 then it follows from the monotonicity principle applied to Z 2 that 
M 2 = for any u- limit point. The monotonicity of Z\ then implies that s — > 1 as r — > 00. 
Lemma 4.2 Consider the dynamical system obtained by restricting the type II system to 
the plane defined by the conditions s = 1 and z = 1. If a solution belongs to the interior of 
the state space for this restricted system then it converges to Pn as r — > 00. 
Proof The restricted dynamical system is identical with that for type II dust solutions. In 



13[ it was proved by using a monotone function derived by Hamiltonian methods that for 

_3_ 

64- 



t — ► 00 the dust solutions satisfy S + — > \ and M 2 — ► Jj. Hence it can be concluded that the 



solution approaches Pn as r — ► 00. 

Lemma 4.3 If a solution lies in the interior of the type II state space then unless it lies on 
the unstable manifold of Pi (and this does occur) the a-limit set consists of the heteroclinic 
cycle described in the statement of Theorem 4.1. 

Proof By Lemma 4.1 we know that any a-limit point satisfies z = 0. Moreover it satisfies 
M 2 = or s = 0. The situation is very similar to that in the massless case treated in pi 
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and the proof may be taken over rather directly. It is only necessary to pay attention to the 
fact that it is the nature of the stationary points in the full massive Bianchi II state space 
which must be taken into account and that the notation is different. 

Suppose that the a-limit set contains a point with s = and M2 7^ 0. Then by Lemma 
4.3 of it contains the endpoint of L\ and either P2 or P3. On the other hand, if it 
contains a point with M2 = then this belongs to the massless Bianchi I state space. Then 
it must contain one of the points Pi, P2, P3, P4, P5 or the endpoint of L\. To prove the 
lemma we may assume that the solution does not lie on the unstable manifold of P\. If Pi 
nevertheless belonged to the a-limit set then this set would have to include points belonging 
to the unstable manifold of Pi other than Pi itself. But these satisfy neither s = nor 
M2 = and so this gives a contradiction. Thus under the given assumptions the a-limit 
set does not contain Pi. If the a-limit set contained a point with M2 = 0, ]£+] < 1 and 
< s < 1 it would contain Pi (in its role as w-limit set for Bianchi type I solutions), leading 
once more to a contradiction. If the a-limit set contained P2 then by the results of the 
appendix it would contain Pi, which is also not possible. Applying Lemma 4.3 of [10] again 
allows points with M2 7^ which are not on the vacuum boundary to be excluded from the 
a-limit set. The straight lines joining P2 to P3 and the endpoint of L\ are excluded as well. 
The conclusion is that the a-limit set is contained in the heteroclinic cycle mentioned in the 
statement of Theorem 4.1. It remains to show that it is the whole heteroclinic cycle. This is 
straightforward to do using the results of the appendix. 

Lemma 4.4 If a solution lies in the interior of the type II state space then it converges to 
the point P\\ as r — > 00. 

Proof Consider any w-limit point with z 7^ 1. Then by Lemma 4.1 this point satisfies s = 1 
and M2 = 0. Any nearby w-limit points must also satisfy these conditions. If any of these 
limit points satisfied z = then P4 and P5 would be w-limit points of the given solution. 
Using the saddle point properties of these points then shows that Pj and Pg belong to the 
w-limit set. Repeating the same argument shows that the endpoint of L\ with s = 1 is an 
w-limit point. The fact that this point is a transversely hyperbolic saddle implies that its 
unstable manifold in the hyperplane s = 1 is contained in the w-limit set. By Lemma 4.2 the 
cj-limit set also contains P\\. Since Pq is a hyperbolic sink this contradicts the assumption 
z 7^ 1. Thus we conclude that the entire u limit set is contained in the plane defined by 
the equations s = 1 and z = 1. The argument just given rules out the possibility of w-limit 
points with M2 = 0. Applying Lemma 4.2 once more shows that the only possible w-limit 
point which does not lie on the vacuum boundary is P\\. Finally the fact that P7 and Ps are 
hyperbolic saddles can be used to rule out points of the vacuum boundary, thus completing 
the proof. 
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5 Type III models 



The physical state space, G, of the LRS type III models is given by the region in R 4 defined 
by the inequalities M 3 > 0, < s < 1, < z < 1, and 1 - - M 3 > 0. 

The coordinates, in terms of s, z, M3), of the various stationary points are the follow- 
ing: (0, So , 0,0),(|, 0,0,0), (1, 0, 0, 0),(1, 1, 0, 0),(-l, 1, 0, 0),(1, 0, 1, 0),(1, 1, 1, 0),(-l, 1, 1, 0), 
(5,0, 0, |),(|, 0, 1, |), where so the same particular constant value of s that appeared in the 
previous type I section. These points are called Pi, ...-P10 (note that they are numbered dif- 
ferently than in |fL0"f , in the massless case). In addition there exist three lines of equilibrium 
points, (—1, 0, K, 0, 0), (0, F, 1, 0, 0), , 1, zq, |), denoted by L\, L2, £3, where K, F and zq are 
constants. The first eight stationary points and the first two lines correspond to points and 
lines of the same name in the Bianchi I system and their coordinates are obtained by append- 
ing a zero to those of the Bianchi I points. The points Pi, P2, P4, Pq, P7, Ps, P9 are hyperbolic 
saddles while P5 and P10 are degenerate, with one zero eigenvalue each. The point P3 is a 
hyperbolic source. The lines L\ and L2 are transversally hyperbolic saddles while the line L3 
is degenerate with two zero and two negative eigenvalues. 

To prove global results about the global properties of solutions it is useful to note the 
existence of the following bounded monotone function 

M 3 = M 3 (2 - S+)- 2 , 

M 3 = 2M 3 [(l-2£ + ) 2 + ft(P + P + )](2-£ + n 1 . (26) 

Theorem 5.1 If a smooth non- vacuum reflection-symmetric LRS solution of Bianchi type 
III of the Einstein- Vlasov equations for massive particles is represented as a solution of ( |l2"|) 
with M2 = 0, then for r — > 00 it converges to a point of the line L3 with z > 0. For r — ► —00 
there exists 

(i) a one-parameter set of solutions lying on the unstable manifold of Pi and 

(ii) a two-parameter set of solutions lying on the unstable manifold of P2 and 

(iii) all remaining solutions converge to P3. 

In all these solutions the scale factor a is monotone increasing at late times. 

Lemma 5.1 If a solution belongs to the interior of the type III state space any a-limit point 

satisfies M3 = 0. Any cj-limit point satisfies S + = ^, M3 = | and s = 1. 

Proof The continuous function M3 on the state space must have a maximum and since its 

gradient never vanishes this maximum can only be attained at points with M3 = 1 — S 2 ^. 

Computing the derivative of M3 along the curve in the (M3, £+) plane defined by this relation 

shows that the maximum value is | and that it is attained when S + = ^ and M3 = |. Now we 

apply the monotonicity principle. Let S be the part of the Bianchi III state space obtained 

by removing the points with M3 = and those with S + — ^ = M3 — | = 0. This is an 

invariant set for the dynamical system. It will now be shown that M3 is strictly increasing 
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along solutions on this set. If S + 7^ ^ or if Q(R + R + ) ^ then this follows immediately 
from (H). If S + = \ and Q(R + R+) = then 

S + = |(M 3 - |) (27) 

This completes the proof that M3 is strictly increasing on S. The monotonicity principle 
then shows that any point in the a-limit set must be in the complement of S and such that 
such that M3 does not take on its maximum value on S there. Hence M3 = there. It also 
shows that any point in the w-limit set must be in the complement of S and that M3 does 
not take on its minimum value there. Hence in the latter case = \ and M3 = |. It follows 
from this that E + — > ^ as r — > 00 and the equation for s then implies that s — > 1. 
Lemma 5.2 A solution which belongs to the interior of the type III state space converges to 
a point of the line L3 with z > as r — > 00. 

Proof Because of the result of Lemma 5.1 it only remains to prove that z tends to a positive 
limit as r — > 00. Note first that the evolution equation for s implies an equation of the form 
(d/dr)(l — s) = (1 — s).F where F = — 6s$] + . As r tends to infinity i 7 — > —3 and a simple 
comparison argument proves that 1 — s(t) = 0(e^ 3+e - >T ) as r — > 00. In particular, 1 — s 
decays exponentially to zero at late times. The evolution equations imply that satisfies the 
equation: 

Si/n = (£+ - |)[(3 - + |(1 + R)] - 2S + (i? + + R) - (M 3 - |)(1 + R) (28) 

Note that R + + R> so that the second term on the right hand side is negative. However it 
is exponentially small at late times since it contains a factor (1 — s) when expressed in terms 
of the matter quantities. In particular Cl/Q — > as r — > 00 and = 0(e eT ) for any e > 0. 
This means that f2 converges to zero slower than any exponential. In other words, VLe tT tends 
to infinity for any e > 0. Suppose that S + > \ for some solution at some time. Then M3 < | 
and the first and third terms in the expression for Q/Q are positive at late times while the 
second term is negative. It will now be shown that the third term decays slower than any 
exponential and thus must eventually dominate the second term. For 

n = (| - Z%) + (| - M 3 ) < (| - M 3 ) (29) 

It follows that &/Q > at late times as long as S + > |. Since it follows from Lemma 5.3 
that Q — > as r — > 00 it follows that for any time tq for which S + (tq) > \ there exists a 
time t > tq with S + = |. When S + = ^ then 

t + = (M 3 - f )[f (1 - it!) - \{R + R + )] (30) 

Now it follows from the evolution equation for z that 1 — z cannot approach zero faster than, 
for instance, e~ T and the same is then true of 1 — R. It can be concluded that the first term in 
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the square bracket on the right hand side of (|3^) dominates the second at late times. Hence 
i — S + must be negative at late times, which in turns implies that z is increasing and that 
it must tend to a positive limit. 

Lemma 5.3 If a solution lies in the interior of the type III state space then unless it lies on 
the unstable manifold of Pi or P2 (and both of these cases occur) the a-limit set consists of 
the point P3. 

Proof Note first that it can be concluded as in the proof of Lemma 3.4 that any a-limit point 
satisfies z = 0. Thus, applying Lemma 5.1, it can be identified with a point of the state space 
for massless type I solutions. Now it is possible to proceed further following the method of 
proof of Lemma 3.4. Consider the boundary of the state space for massless type I solutions. 
The point P3, being a hyperbolic source in the type III state space, can be excluded as an 
a-limit point of a solution of type III. It is then possible to successively exclude points of 
the boundary as in the proof of Lemma 3.4. The facts which need to be used are that all 
a-limit points satisfy M3 = and z = 1 and that the points P4, P5 and the endpoint of 
L\ are a hyperbolic saddle, a non-hyperbolic saddle topologically equivalent to a hyperbolic 
one and a transversely hyperbolic saddle, respectively. At this stage it can be concluded 
that all a-limit points of solutions of type III are either Pi, P2 or points of the unstable 
manifold of Pi. For all other points of the interior of the massless type I state space lie on 
solutions which converge to the hyperbolic source P3 in the past time direction, and so are 
excluded. It remains to examine what happens in a neighbourhood of the points Pi and P2, 
which are both hyperbolic saddles. The unstable manifold of P2 in the type III state space 
is three-dimensional and so there are solutions which converge to P2 as r — > — 00. Any other 
type III solutions which had P2 as an a-limit point would have to have a-limit points on the 
stable manifold of P2, which has already been excluded. Hence solutions of type III which do 
not converge in the past to P2 cannot have P2 or a point of its unstable manifold as a-limit 
points. Thus the only remaining possibility is that solutions lie on the unstable manifold of 
Pi and converge to that point in the past. Since the unstable manifold is two-dimensional, 
solutions of this kind exist. 

The results of Lemma 5.2 and Lemma 5.3 together imply all the results of Theorem 5.1 except 
the last directly. The statement about the scale factor a follows from the fact, derived in the 
course of the proof of Lemma 5.2, that £+ < i at late times. 

6 Concluding remarks 

In this paper we studied the dynamics of solutions of the Einstein- Vlasov equations which are 
locally rotationally symmetric, reflection-symmetric and of Bianchi types I, II and III. The 
initial singularities are of four types. There are isotropic singularities which, in the dynamical 
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systems description used in this paper, are those which converge to the point Pi as r — > — oo. 
The general theory of isotropic singularities developed by Anguige and Tod pj implies as a 
very special case the occurrence of isotropic singularities in Bianchi models with collisionless 
matter and information about how many there are. They only developed the theory for 
massless particles and so in order to apply to the situations considered here it would have 
to be generalized to the massive case. There are barrel singularities which occur in types I 
and III but not in type II. In the dynamical systems picture these are the solutions which 
converge to P2 as r — » —00. Fluid models with corresponding symmetries never have barrel 
singularities and so this is a peculiarity of collisionless matter, both in the case of massive 
particles studied here and that of massless particles studied in [10|. There is the generic case 
in types I and III, which concerns solutions which develop from an open dense set of initial 
data for each of these Bianchi types. These solutions have a cigar singularity and converge to 
-P3 as r — ► — 00. Finally, there are the generic solutions of type II, which have an oscillatory 
initial singularity. 

As far as the late time behaviour is concerned, it is tempting to speculate that behaving 
like a dust model at late times in an expanding phase may be a general feature of solutions of 
the Einstein- Vlasov equations with massive particles. We know of no counterexample to this. 
For the solutions of types I and II treated in this paper it has been proved to be true. For type 
III the situation appears to be delicate and the occurrence of degenerate stationary points 
of the dynamical system may require an application of centre manifold theory in order to 
determine details of the asymptotics. A possible criterion for detecting cases where there may 
be trouble is as follows. Consider a dust solution which is a candidate for the asymptotic state 
of solutions of the Einstein- Vlasov equations. If each eigenvalue of the second fundamental 
form of the homogeneous hyper surf aces, when divided by the mean curvature, is bounded 
below by a positive constant in the dust solution then it is a strong candidate. Otherwise 
difficulties are to be expected. This criterion gives a positive recommendation for types I and 
II and a warning for type III. Thus at least for the models investigated in this paper it is a 
good guide. Using the information on dust models in chapter 6 of jDJ it also gives a positive 
recommendation for types I, II and VLq without the need to restrict to the LRS case. 

In this paper a dynamical system has been set up for all LRS Bianchi models of class 
A as well as for Kantowski-Sachs models and the type III models, which are of class B. We 
expect that techniques similar to those used here can be applied to analyse Kantowski-Sachs 
models and LRS models of type VIII and IX. An important feature of all these LRS models 
is that the Vlasov equation can be solved exactly. This is also true of general Bianchi type 
I models. Some limited results on the dynamics of Bianchi type I solutions of the Einstein- 
Vlasov equations which are reflection-symmetric but not necessarily LRS were proved in 

A heuristic analysis of reflection-symmetric type I models was given using Hamiltonian 
techniques in 0, where there are also interesting remarks on the general Bianchi I case. It 
would be very desirable to have a mathematically rigorous implementation of the ideas of M . 
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What can be done in cases where the Vlasov equation cannot be solved exactly? If, as 
already speculated above, the late time evolution resembles that of a dust solution and if the 
dust solutions are asymptotically LRS then it may be possible to give a good approximation 
to the solution of the Vlasov equation in that regime. There is one drawback of this idea 
as a general tool for Bianchi class A models. Unfortunately there are no LRS spacetimes of 
Bianchi type VL> In the case of fluids there exists a special class of Bianchi VLj spacetimes 
which is often characterized by the rather opaque statement that = 0. These spacetimes 
do have a simple geometric characterization which will now be explained. Every Bianchi 
class A spacetime has a discrete group of isometries whose generators simultaneously reverse 
two of the invariant one- forms on the group. The special class of Bianchi VIo solutions 
can be characterized by the existence of an additional isometry which reverses just one of 
the one-forms. It is possible to consider solutions of the Einstein- Vlasov equations with the 
corresponding type of symmetry. We are, however, not aware that the Vlasov equation can 
be solved exactly in these special spacetimes. If it could then this might fill the apparent gap 
in the strategy just suggested. 

The oscillatory behaviour observed near the singularity in type II models appears at first 
sight to indicate that collisionless matter does not fit into the analysis of general spacetime 
singularities by Belinskii, Khalatnikov and Lifshitz Q]. On the other hand, the fact that 
in the analysis of Misner |?]] using a time-dependent potential we see the phenomenon of 
walls moving too fast to be caught suggests that the oscillations might go away in general 
models. This issue requires further work. It could turn out that collisonless matter generically 
becomes negligible near the singularity, as originally stated for fluids in 0]. 

To conclude, we mention some further interesting open problems. What happens in the 
case of a model with two species of particles, one massive and one massless? Of course this 
could be thought of as a simple cosmological model incorporating both baryonic matter and 
the microwave background photons. It is related to the two-fluid models which have been 
analysed in ||. Mixtures of fluids and kinetic theory could also be considered. We have seen 
that the Einstein equations with collisionless matter as source may behave very differently 
from the Einstein equations with a fluid source at early times (and also at late times in the 
massless case) . Under what circumstances are there intermediate stages of the evolution with 
collisionless matter which can be well described by a fluid? Since the point Pi is a saddle there 
are obviously solutions which approach this point and then go away again but is there more 
that can be said about this issue? What can be said about inhomogeneous models? In || 
Rein analysed the behaviour at early times of solutions of the Einstein- Vlasov equations with 
spherical, plane and hyperbolic symmetry and massive particles. He identified open subsets of 
initial data for these symmetry types with a singularity resembling the generic LRS solutions 
of types I and III. There is an overlap between the results of Q and those of the present 
paper. It could be illuminating to attempt a common generalization of these. In any case, it 
is clear that one of the central challenges of the future in the study of cosmological solutions 
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of the Einstein- Vlasov equations, or indeed the Einstein equations coupled to any type of 
matter fields, is to develop techniques which apply to inhomogeneous problems. A thorough 
understanding of the homogeneous case is likely to be an invaluable guide in addressing it. 



A Appendix 

In this appendix some general procedures which are useful in determining limit sets of solu- 
tions of dynamical systems will be outlined. Let 7 be an orbit of a dynamical system and 
p a stationary point. We will discuss only cj-limit sets, but corresponding statements about 
a- limit sets follow immediately by reversing the direction of time. We consider the following 
three statements which may or may not be true for given choices of 7 and p. 

1. p is an cj-limit point of 7 

2. 7 lies on the stable manifold of p 

3. there are w-limit points of 7 different from p which are arbitrarily close to p and lie on 
the unstable manifold of p 

4. there are w-limit points of 7 different from p which are arbitrarily close to p and lie on 
the stable manifold of p 

In the body of the paper we frequently use certain relations among the statements above 
which hold under various assumptions on the nature of the stationary point p. Whatever the 
stationary point, it is always true that the statement 1. is implied by any of the statements 
2., 3. or 4. This is a consequence of the elementary fact that the w-limit set is closed. Now 
suppose that p is a hyperbolic stationary point. In this case, if 1. is true and 2. is false then 



both 3. and 4. are true. This follows from Lemma Al of |1C]. Combining these statements 
we see that for a hyperbolic stationary point there are two mutually exclusive cases under 
which 1. can hold. Either 7 lies on the unstable manifold of p or the w-limit set contains 
points of both the stable and unstable manifolds of p arbitrarily close to p. In particular, 
if p is a hyperbolic source then it cannot be in the w-limit set of 7 and if p is a hyperbolic 
sink and p is in the w-limit set of 7 it is the whole w-limit set. If we already have some a 
priori information about where w-limit points can lie (due, for instance, to the existence of 
a monotone function) then this gives more information about where the points on the stable 
and unstable manifolds whose existence is guaranteed by the general statements above can 
lie. 

Next we consider the case of transversally hyperbolic stationary points. Suppose that p 
belongs to a manifold of stationary points of dimension d. (Only the case d = 1 occurs in 
this paper.) These points have a zero eigenvalue of multiplicity d. If all other eigenvalues 
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have non-vanishing real parts then the stationary point is called transversally hyperbolic. 
(Depending on the signs of the eigenvalues the manifold of stationary points is called a 
transversally hyperbolic source, sink or saddle.) By the reduction theorem (Theorem Al of 



[10 1) each of these points lies on an invariant manifold and the restriction of the flow to each 
invariant manifold is topologically equivalent to that near a hyperbolic stationary point. The 
arguments for a hyperbolic stationary point adapt easily to give analogous statements for 
transversally hyperbolic stationary points. In applying these results we can essentially ignore 
the directions along the manifold of stationary points. 

Finally we consider certain other non-hyperbolic stationary points. A result of the type we 



need was proved in Lemma A2 of [1C] but we would like to formulate the statement in a more 



transparent way here. Consider an isolated stationary point p with a trivial stable manifold 
and a one-dimensional centre manifold. Using the reduction theorem we see that the unstable 
manifold divides a neighbourhood of p into two parts on each of which the restriction of the 
dynamical system is topologically equivalent to the restriction of a dynamical system with a 
hyperbolic stationary point. Whether the latter system has a saddle or a source depends on a 
certain sign condition. This condition may be different for the two halves. In the dynamical 
systems considered in this paper the only example of this is provided by the point P5. Only 
one of the halves belongs to the physical part of the state space and in that half the sign is 
such that a saddle is obtained. The result of these considerations is that for the arguments 
in this paper P5 may be treated just as if it had been a hyperbolic saddle, with the centre 
manifold taking over the role of the trivial stable manifold. 
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